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A slave rotor-Hartree Fock formalism is presented for studying the properties of the p-d model 
describing perovskite transition metal oxides, and a flexible and efficient numerical formalism is 
developed for its solution. The methodology is shown to yield, within an unified formulation, the 
significant aspects of the rare earth nickelate phase diagram, including the paramagnetic metal state 
observed for the LaNiOs and the correct ground-state magnetic order of insulating compounds. It is 
then used to elucidate ground state changes occurring as morphology is varied from bulk to strained 
and un-strained thin-film form. For ultrathin films, epitaxial strain and charge-transfer to the apical 
out-of-plane oxygen sites are shown to have significant impact on the phase diagram. 
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Understanding the unusual electronic behavior of 
transition-metal oxides has been a long-standing ques- 
tion in condensed matter physics and interest has in- 
tensified following the demonstration^ that the materials 
could be used as components of atomically precise ox- 
ide heterostructuresi^"— The theoretical challenge posed 
by the materials is to treat simultaneously the strong 
local correlations in the transition metal d-orbitals and 
their substantial hybridization with oxygen p orbitals. 
In some systems, the p orbitals can be integrated out 
and the physics represented in terms of a (possibly multi 
band) Hubbard model representing the d orbitals only, 
for which many theoretical methods are available^. How- 
ever, in many cases the charge transfer between p and d 
orbitals is large enough that the p states cannot be ne- 
glected. This "negative charge transfer energy" regime^ 
has been less extensively studied. While much useful 
information has been provided by density functional the- 
ory (DFT)1 and extensions such as DFT+U^— , hybrid 
hmctionalsi^i^, and DFT+DMFTiii^, these methods 
are computationally intensive, so that the large super- 
cells required for long-period ordered phases are difficult 
to study. Furthermore, the variety of experimental bulk 
and superlattice configurations and of many-body phe- 
nomena emphasizes the need for a model-system treat- 
ment that encapsulates the essential physics so the im- 
portance of different contributions can be disentangled. 

The rare earth nickelatesi^, ReNiOs, are an impor- 
tant case in point. According to standard valence count- 
ing arguments, Ni is in the 3-1- oxidation state with 
a d'' valence configuration. However, photoemission 
experiments2i reported that the dominant GS configu- 
ration is d^L, in agreement with unrestricted Hartree- 
Fock (UHF)2ai2i j^jj^ DFT-hUl^-ii, and DFT+DMFTil 
approaches. Additionally, O Is x-ray absorption exper- 
iments have found significant hole concentration on the 
oxygen p orbitals^ir— . This large degree of charge trans- 
fer places the material near the negative charge-transfer 
regime^i^. As Re is varied across the lanathanum row 
of the periodic table, the ground state (GS) of bulk ma- 
terials changes from paramagnetic metal (PM-M) to cor- 
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FIG. 1: (Top to bottom) Bulk, 2-, and 1-layer phase diagram 
on the u-Sdpd plane. Slight charge ordering accompanies the 
spin ordering for Sdpd ^0. lt/t/0 denote large moment/small 
momont/non-magnetic sites of the long-period magnetic or- 
dering patterns detailed in the text. In the lower two panels 
the light dotted line indicates the phase boundary obtained 
without the energy shift on the apical oxygen. 



related insulatori^i^. The correlated insulator phases 
exhibit a rock-salt-pattern lattice distortion^"— and a 
nontrivial long-period magnetic ordering^^"— . In ultra- 
thin films a metal-insulator transition which is appar- 
ently unaccompanied by rock-salt-pattern lattice distor- 
tion occurs as film thickness and strain are varied^"—. 
Understanding how these apparently different transitions 
can occur within a single formulation is an important 
open theoretical challenge^"— i2£ii2ri^. 

In this paper, we develop a mean-field (MF) approach 
based on a combination of slave-rotor (SR)^"— and 
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Hartree-Fock (HF) methods, along with an efficient and 
flexible numerical formalism for solving the MF equa- 
tions. The method is powerful enough to permit the ex- 
amination of the large supercells needed to investigate 
long-period ordering patterns in the context of realistic 
crystal structures. The results reconcile the bulk and 
film phase diagrams and allows us to identify the key 
role played by the oxygen degrees of freedom and the 
lattice distortions. 

Figure [1] summarizes our key new results. It presents 
phase diagrams in the space of Ni charging energy U 
and amplitude Sdpd of rocksalt-type lattice distortion 
(defined more precisely below). The upper panel shows 
that, in the absence of oxygen breathing distortion, bulk 
materials are PM-M at any value of U, consistent with 
experiments^. As the lattice is distorted, transitions oc- 
cur, first to a metallic magnetic phase and then to a 
magnetic insulator. The nontrivial ordering wavevector 
found in experiment^"— is correctly obtained as a 3D 
16-formula-unit magnetic pattern of the type tfOJJ-O, de- 
noting the z-projcction of magnetic moment detailed be- 
low. The 4 fl" and 4 JJ- Ni sites are surrounded by 6 spin-0 
Ni sites, which in our interpretation form singlets with its 
6 oxygen neighbors. Each of the 8 spin-0 Ni sites are con- 
nected to 3 11 and 3 ^ Ni sites via oxygen bonds. States 
with the ttii pattern are unstable against PM or -ff 4 
solutions, and indeed are not observed in experiments^. 
In agreement with symmetry-based arguments^, slight 
charge ordering always accompanies the spin ordering. 
For reasonable parameter values the change in electronic 
structure across the Re series is properly accounted for, 
with one exception: the magnetic metallic phase found in 
theory is not observed in experiment. This is discussed in 
more detail below. The lower panels show the evolution 
of the phase diagram with film thickness M and (bottom 
panel) with applied strain, revealing that in ultra-thin 
films an insulating phase can occur even in the absence 
of rocksalt-type lattice distortions. 

Our theoretical approach begins from a p-d lattice 

model of the form H ^ Y.k + J2i i^i + Ji), with 
k and / sum over crystal momenta and Ni sites, respec- 
tively. The parameters of this model may be obtained 
from, e.g. maximially localized Wannier fits to DFT 
results^, but the precise form is not important for this 
paper (see also Ref. [10 for another example of the insen- 
sitivity of results to the precise p-d model band parame- 
ters). We adopt a simplified scenario in which we retain 
the Ni Cg and O 2p(7 orbitals, with nearest-neighbor p-d 
and p-p hopping. The difference between the bare Ni-d 
and 0-2p energies Ed-Sp is important, as discussed below. 

The structure of the ReNiOa materials is derived from 
the ideal ABO3 cubic perovskite structure which is a lat- 
tice of corner-sharing oxygen octahedra, each contain- 
ing a Ni site at the center. The crystal structure of the 
actual materials is distorted from this structure by ro- 
tations of the octahedra which are not important for 
our purposes and, in the insulating cases, by a two- 
sublattice distortion in which adjacent Ni's have signif- 



icant different Ni-0 bond lengtha^^"— . To incorporate 
the bond disproportionation in the Hamiltonian under 
the cubic approximation, we scale the hopping accord- 
ing to the Harrison rulo^, tpd = tpd (j- + ^'^pd/'Ppdj 
and tpp = t^pp (1 + 6dpp/dpp) ^, with d^^ = 1.95A and 
dpp — y/2dpj^. An additional effect may occur in layered 
structures. Liu et al.^ showed that the presence of Al 
at the interface would deplete holes on the out-of-layer 
oxygen sites linking Al and Ni, raising the charge-transfer 
energy from those apical sites by '--^leV. We model the M- 
layer 2D structures using supercells with Al NiOs units 
in the z-direction, terminated on both ends with apical 
oxygen sites whose Cp are shifted by -leV. For reference 
we also performed some calculations without such shift. 

The interaction terms Ui and J; contain the on-site 
repulsion and the Hund's interactions in the rotationally- 
invariant Slater-Kanamori formi. 

Ul = U ^ rimcrnm'cr' (1) 

(m(T)>(m'cr') 



cr 771 

- j^(^^naanb^a+dl^dl_^dbada-a^ 

+ J {dlAi'^bidt^ + ^■'^■) (2) 

In this expression, J/ differentiates only states with the 
same occupancy but inequivalent configurations. 

To treat the on-site Coulomb interactions, we adopt 
the slave rotor (SR) approach^"—. The method was orig- 
inally applied to the Hubbard model and then to other 
d-only models with Hund's-like interactions^!^ under 
the j << U approximation. We apply it here to the 
p-d model, noting that for RNiOs j^leV is much smaller 
than either the d-d repulsion and the electron bandwidth. 
For each Cg site, the approach introduces an auxiliary SR 
field, 6i E [0, 27r), and decomposes electron operators as 

'^Zqct ^ fiaa^^^^- "^^^^ consistency of SR state and d- 
occupancy is enforced by the constrain 




The spectrum of C// ' cc Lf. We treat the large Ui with 
9 and the small j-scale J; with using the weak cou- 
pling Hartree-Fock approximation. That is, we solve 
the p-d model with the singlc-Slatcr-dcterminant ansatz 
\MF) = \p, f)\9). We follow previous SR applications us- 
ing Lagrange multipliers, hi. Along with constraint Eq.[3] 
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and up to a constant, the system of equations reads 

rif) ■ -3, 



I laa 

J2{e-''')V,o.pU,Jic.a + h.c. 

leaa 



(4) 



le6(T 



He 



hiLi 



(5) 



where V^^ is given by the product of t^^; and the Slater- 
Koster orbital symmetry factors in Table I of Ref. [52|. 

We solve the MF equations without any restriction 
other than the single-Slater-determinant assumption. A 
3D Bravais lattice of Nii6 048 supercells is required for 
unrestricted modeling of the (1/2,0,1/2) pattern with 
respect to the ortliorhombic unit cell. To systemat- 
ically capture the effects of dimensionality, we model 
the M-layer structures using a 2D Bravais lattice of 
Ni4MOi2M+4 supercells, which is connected to the 3D 
Nii6048 supercell for M — >• oo. The system is solved 
by the standard T=0 iterative procedure for up to 
65536 supercells. Without much optimization, the worst 
case initial condition with bi-partite charge-, orbital- 
, and magnetic-order would converge into a PM-M or 
(1/2, 0, l/2)-ordered insulator within 30 cpu hours on an 
Opteron-2350 cluster. 

We perform bulk calculations for {vL,ed) such that 
{ugg ) ~ 2 as found in UHF with parameters fitted to pho- 
toemission spectrum2£, DFTi^^ii, and DMFT+DFTJ^. 
For the reference set of ?7 = 5.3, j = 1, tp„,(i^2_„2 = 1-5, 
tpp ~ 0.5, Apd = Cd — Cp + U ~ 0.3, the undistorted 3D 
lattice has (ng^) ~ 1.95, slightly less than those of the 
above calculations. To isolate dimensionality effects, we 
use the bulk's {u,ed) pairs, in addition to the aforemen- 
tioned apical Cp shift, for layered calculations. 

The left panel of Fig. [2] shows representative results for 
the energies of different locally stable phases computed 
at interaction (U,j)=(5.3,l) as a function of {SdpdY ■ We 
see that in agreement with experiment, the undistorted 
lattice is a PM-M. As oxygen breathing disproportiona- 
tion 5dpd is increased a transition to a period-4 metallic 
magnetic state occurs. At yet larger distortion a metal- 
insulator transition occurs. The ferromagnetic states 
found in UHF^i^i , DpT+XJi^-ii and DFT+DMFTii be- 
come stable only at larger 5dpd and gain less energy than 
does the experimentally observed ordering pattern. Also, 
E^—Sdpj^ implies an equilibrium distortion ultimately 
determined by the aharmonic lattice restoration force be- 
yond the scope of this study. 

The dashed and dotted lines in the right panel of Fig. [5] 
are the SR renormalizations of the p-d hopping (Eq. 2]) 
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FIG. 2: (Left) Stable solutions' MF energy. (Right) rotor 
renormalization and insulating gap for the bulk system at 
[/ = 5.3. 



on the two sublattices of the distorted structure. We see 
that the Coulomb renormalizations are not large, and 
are only weakly dependent on sublattice. This renor- 
malization has a moderate effect on physical properties, 
e.g. Fermi velocity renormalization at (U,j) = (5.3,l) is 
''^5.3,1/^^0,0 ~0.65. The relatively modest SR renormaliza- 
tions confirms that near the "negative charge transfer" 
region of the phase diagrani^ the on-site Coulomb corre- 
lations do not drive the phase transition in an orthodox 
way, in particular the physics is far from the (e*^) 
Brinkman-Rice insulator. Even then, the renormalized 
kinetic terms, combined with lattice distortion and the 
magnetic order opens a gap in the density of states. We 
stress that SR renormalization (e*^) due to U is required 
for insulation within physical parameter range. 

Parameters scaled for different Re bulk materials^^"— 
are marked on the phase diagram in top panel of Fig. [T] 
It captures all experimentally observed phases over the 
Re series, with the exception of Re=Pr, which is pre- 
dicted to be a metal instead of an insulator, albeit with 
the correct charge and magnetic order. This can be un- 
derstood by noting that our approach underestimates the 



insulating gap, which is 160meV at ^^=2.5% compared 

to -200meV in a 2-sitc DFT+DMFT studyii, and that 
Re=Pr is also extrapolated to be a metal from those 
Re=Lu results. 

In agreement with experiments, the tfOJJ-O insulator 
(top panel of Fig. [T]) has a (1/2,0,1/2) magnetic structure 
with respect to the orthorhombic unit cell, a rock-salt 
charge ordering pattern, but without orbital ordering. 
Let us identify the Ni sites with longer/shorter Ni-0 bond 

''^=2.5% 



length as Nij^^^^g and Nig-^Qj.|. and consider the ■ 
solution. The occupancy at those sites are n 



long 



2.03 



and ng]-^pj.^=1.90 . This disproportionation can be un- 



derstood by measuring the Ni-0 hybridization order pa- 
rameter 5, = Eeo..yLe-'''pl+eJio.a + h.c. We found 
<°shortl 



that 



1.4, which is greater than expected 
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FIG. 3: Results for u — 7. (Top) d-occupancy for the 1-layer 
structure. (Bottom) raiong magnetic moment for bulk, middle 
layer of 3-layer, 2- and 1-layer structures. 



from the corresponding Stpd/t^^ ^10% hopping modula- 
tion. This shows that Nigj^Qj.^- is strongly hybridized with 
the oxygen bands while Nij^^^-^g retains its d^ characteris- 
tic. The magnetic moments are 771-^^^^=1.1, which agrees 

with experimen t^^'^^ , and "^ghort"'^' charge-transfer 
insulators, the effective physics of strong p-d hybridiza- 
tion is a strong antiferromagnetic Sd ■ Sp spin correlation. 
In this context, having nig^^^^ = in the single-Slater- 
determinant approximation is evident of a d^-high-spin 
state forming a singlet with its oxygen neighbors. The 
formation of p-d singlets has also been reported in a two- 
site DFT-t-DMFT study which predicted FM GSii. 

We now discuss the effects of dimensional confinement. 
The bulk and M > 3 layers structures have similar phase 
diagrams, with increasing ordered and insulating region 
for decreasing M and increasing U. In particular, the 
lower panel of Fig. [3] shows that PM-M exists for M > 3. 
The phase diagram of 2-, and 1-layer ultra-thin films are 
shown in the lower panels of Fig [1] We see that, for 
M < 2 layers, the ground state is always magnetically 
ordered, even in the absence of the bond disproportion- 
ation. Further, the distortion amplitude required for the 
transition to the insulating phase decreases with increas- 
ing u and decreasing M. At (u,j)=(5.3,l), the metal- 
insulator boundary (^^)a//t is reduced from 1.9% to 

0.96% and 0.44% for 2- and 1-layer, respectively. A com- 
parison between the solid and dotted lines in Fig. [T] shows 
that charge transfer to the out-of-plane orbitals^ greatly 
enlarge the insulating regions. 

To study the effects of epitaxial strain, we focus on 
1-layer structure for concreteness. We include the strain 



effects by scaling hopping integrals according to the Har- 
rison rule^ for the geometry specified in Ref. [sol . which 
introduced compressive (tensile) strain with LaSrA104 
(SrTiOs) such that the lattice parameters are a=b=3.769 
(3.853)A and c = 3.853 (3.790)A . Breathing-mode dis- 
tortion is then introduced as discussed above. In agree- 
ment with experiments on ultra-thin layer o'^^'^^i'^^ . Fig. [1] 
shows that compressive (tensile) epitaxial strain enlarges 
(shrinks) the insulating region. Under realistic values 
of U and compressive stain, ordered insulator is possi- 
ble for Sdpd=0, agreeing with experiments which found 
breathing-mode distortion only under tensile strain^. 
Note that further shift in apical oxygen Cp also decreases 
Sdpd required for insulation. 

Figure [3] demonstrates an orbital polarization of ~5%. 
We also found that the polarization decreases (increases) 
with compressive (tensile) strain, in agreement with a 
DFT-fU studjiii^. For example, LaSrA104 (SrTiOg) 
substrate changes the polarization to ~4 (6)%. While 
DFT+U predicted FM order, our M=l,2 GS have an 
AFM in-plane order of fft-IJ-ii which differs slightly from 
the bulk's ffOJJ-O pattern by having small nishort^Q which 
decreases with increasing distortion. The 2-layer struc- 
ture has an additional magnetic transition between FM 
and AFM ordering in the z-direction. 

Lastly, we compare our results with other theoreti- 
cal works. In our formulation, the bulk PM-M devel- 



ops into lyO-O-O-I for 



>1.9% compared to UHF's 



>7.5% transition from ferromagnetic to MCO state 

with similar pattern^. The pattern is also similar to 
the j/U > 1 S-SDW state in the d-only model^, but we 
have j/U < 1/5 and a bond and charge disproportiona- 
tion with d^^^LcP^^L involving strong pd hybridization. 
Unlike the d-only model, our results are sensitive to di- 
mensional confinement. Even though different DFT im- 
plementations have been employed to capture the bulk 
LaNiOa PM-Mi^ and LaNiOa/LaAlOa ordered insulat- 
ing layer—, the demonstration of all bulk and layered 
phases by a single formulation had been elusive to the 
best of our knowledge. 

Our formulation bridges the gap between HF-like ap- 
proaches and the expensive cluster DMFT. The results 
connect all ReNiOs delicate phases in bulk and layer 
form and also provide detailed insights about them. This 
demonstrates a viable pathway of treating systems with 
different superstructures as well as compounds such as 
Fe and Co oxides with strong p-d hybridization and large 
number of partially filled strongly correlated orbitals. For 
example, Sr2Fe04 has also been shown to exhibit strong 
hybridization and non-trivial magnetic order—. 

We thank G. A. Sawatzky, H. Chen, H. T. Dang, R. 
Fernandez, S. Park, and D. Zgid for helpful discussions. 
This effort is supported by US National Science Founda- 
tion under grant NSF-DMR-1006282. 
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